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CHAPTER 1 


INTRODUCTION 

Given a model potential, it is possible, atleast in 
principle, to obtain an equation of state from statistical 
thermodynamics and for this purpose many model potentials 
have been proposed starting with simple hard sphere potentials. 

Several of these model potentials, like the hard sphere 
and Lennard -Jones potential, have been extensively studied, 
Virial coefficients and transport properties are evaluated 
directly from the potentials and thermodynamic properties 
through an equation of state. However, the hard sphere model 
and the L.J. potential give reasonably accurate results only 
for simple monatomic substances. But for polyatomic molecules 
they fail to give satisfactory results*, some of the reasons 
for which are discussed in the latter chapters. Since we are 
concerned with polyatomic molecules in this thesis we are 
looking for a more representative model potential. The m-6-8 
potential proposed by Hanley and Coworkers is one of the 
obvious choices for this purpose. A detailed discussion of 
this potential is presented in the next chapter. 

Even if we have a good model potential, it would be of 
little consequence if an equation of state can not be obtained, 
from that model. Most often this is the crux of the problem 
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and until this is achieved one does not really know which, is 
a good model potential. Various thesis based on perturbation 
technique and variational methods are currently being tried out. 
The so-called computer experiments played a significant role in 
the development and assessment of the accuracy of these thesis. 

Though the results of perturbation calculations have been 
satisfactory for a large number of potentials in a given tempe- 
rature range, the mathematics involved in the expansion of the 
radial distribution function in powers of perturbation potentials 
are cumbersome. Often this series has to be terminated with the 
second order terms and the errors involved in this is significanl 
at low temperatures or when attractive forces are important. Bui 
higher order perturbation theories are exceedingly complicated. 

Hence in this work the variational method (discussed in dets 
in Chapter 4) is used to obtain an equation of state for the 
m -6~8 potential. Using the potential parameters applicable to 
pure CO 2 the equation of state and the vapour liquid coexistence 
curve have been obtained and compared with experimental measured 
values. It is found that the agreement is quite satisfactory. 
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CHAPTER 2 
m-6-8 POTENTIAL 

2.1 Introduction 

H.J.M. Hanley and MarKlein [1] of MBS, Washington, 
proposed a four parameter potential function based on the 
semi-theoretical m-6-8 potential. 

The semi-theoretical potential is 

0 (r > %rr! - -f t 2 - 1 - 1 * 

r r 

where A, m, Cg and Cg are some constants. 

The model proposed by Hanley and Klein is 

/>(*■) = e [f±| v (*/r) m d m - a^I|=8l (a/r) 6 d 6 

- v (cr/r) 8 d 8 ] (2.1.2) 

In reduced form (2) becomes 


AO . ,« 


- V (l/r* ) 8 (2 .1.3) 


In equations (2.1,1), (2.1.2) and (2.1.3) 
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(i) a is defined by the condition 0(cr) = 0 

(ii) d is the ratio of r and a 

(iii) r m is defined by 0* (r ) = © 

(iv) e is given by <p(r m ) = -e 

(v) 0* = 0/e 

and 

(vi) r* = r/r m = r/d a 

The four parameters in this model are m, v and a (or Y m ). 

1 m' represents the strength of the repulsion while v represents 
the strength of the inverse eighth power attraction between the 
molecules. In the following we will refer to (2.1.2) and (2.1,3) 
as the m-6-8 potential for simplicity. 

2,2 Selection of P arameters for m-6-8 Mode l 

The four parameters m, v , a and e/k are calculated either 
by fitting experimental viscosity data or from collision inte- 
gral calculations. The former is simple but not accurate while 
the latter is comparatively more accurate but laborious. 

Ely and Homley [2] used the two procedures simultaneously to 
optimise the time and accuracy. The actual details are described 
in ref. [2], 

Normally v is selected for each m such that m- v (m-8)> 0 

while m is selected between 9 and 18. The following table 
provides parameters of the m-6-8 potential obtained by Ely and 
Hanley for some common gases. 
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TABLE 1 


Gas 

m 


Or(S) 

e/k (°K) 

Argon 

ii 

3 

3.292 

153.0 

Krypton 

ii 

3 

3.509 

216.0 

Xenon 

ii 

3 

3.841 

295.0 

Oxygen 

10 

1.0 

3.450 

109.5 

Methane 

11 

3.0 

3.680 

168.0 

Nitrogen 

12 

2.0 

3.540 

118.0 

Garb ond ioxid e 

14 

1.0 

3.680 

282.0 


2.3 Applications of m-6-8 Potentia l an d Results 

This model potential had been applied to Argon, Krypton, 
Xenon, Nitrogen, Methane and Carbondioxide to check its 
'utility and validity* [3]. Properties like viscosity coefficien 
self diffusion coefficient, thermal conductivity, virial coeffi- 
cients and the isotopic thermal diffusion factor were calculated 
and the results were found to be in excellent agreement with 
experimental values over a wide temperature range [3]. 

Caligaris and Rodriguez [4] calculated the reduced third 
virial coefficient for Argon assuming pairwise additivity 
with m = 11, a = 3.292 j? # r ffi = 3.667 5 and e/k = 153.0 °K. 

They compared the result of this calculation to the experimental 
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values of Argon and found a wide deviation between them. However, 
Hanley and McCarty [5] corrected their results using non-additive 
terms and found a good agreement with experimental results. 

Thus the spherically symmetric m-6-8 model proves to be 
reasonable for non-spherical polyatomic molecules like Carbondi- 
oxide and Methane also. But to account for the orientation 
dependence in such molecules, angle dependent terms should be 
added as shown latter. Self diffusion coefficients of Methane were 
calculated using the molecular dynamic technique [6] for m-6-8 
potential and agreement with the experimental data is within 
.estimated errors. Collision integrals and dielectric second 
virial coefficients were evaluated for a non-spherical m-6-8 
potential function for Carbondioxide [7]. Agreement between 
calculated and experimental values was 'generally satisfactory' . 

These applications and their results show the capability of 
m_6-.8 potential functions in the prediction of physical and trans- 
port properties for spherical as well as non-spherical molecules 
and encourage us to apply this model for the prediction of phase 
equilibrium properties. 

2.4 Superiority of the m-6-8 Model 

The three parameter potentials such as Lennard-Jones 
potential may be compared to the m-6-8 model. The three 
parameter potential functions have three major limitations : 



7 


(i) A set of parameters selected to fit a property 

in a given temperature range cannot be generally relied upon 
to predict the property in some other temperature range. 

(ii) A set of parameters which predict one kind of property 
do not predict another kind of property correctly. 

(iii) The parameters do not always agree with values 
obtained from independent information such as scattering 
experiments. 

On the other hand the m-6-8 potential has been used to 
correlate a given property over a wide temperature range and 
is able to predict properties of different types with a single 
set of parameters [3]. Thus the addition of the inverse eighth 
term to a three parameter function substantially improves the 
relation between theory and experiment and also with this addition, 
the limitations mentioned earlier are essentially removed. It is 
worth noting that the repulsive exponents in this model came out 
to be very close to the results of scattering experiments. 

2.5 Weaknesses of m-6-8 Model 

Th ough m-6-8 potential is an effective tool to handle 
simple monatomic gases, it is not free of limitations. The 
two important limitations of m-6-8 model are discussed below. 
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(i) Dispersion coefficient s 

Prom (2.1.1) and (2.1.2) the dispersion coefficients Cg 
and Cg are given by 

°6 = 5T<j( m -v C m_8 )) aS (2.5.1) 

°8 = ^ (2.5.2) 

Prom (2.5.1) and (2.5.2), it is evident that dispersion coeffi- 
cients are strong functions of a; a is obtained from experiments 
and since a is raised to sixth or eighth power, it can introduce 
considerable uncertainty in Cg and Cg. 

To reduce this uncertainty, one immediate suggestion is 
that more parameters should be added to the model. This view 
can be considered by referring to a potential function due to 
Barker [8] r in which he has added additional parameters to the 
repulsive term. 

The m-6-8 potential with m = 11 and v = 2 is, however, 
very close to Barker's multiparameter potential for which Cg 
is fixed at the 1 correct* value. It would thus appear that 
adding more parameters to m-6~8 potential would not necessarily 
remove the uncertainty. Instead of keeping Cg and Cg as con- 
stants for all the temperature and density conditions we can 
change them by varying the parameter a to fit the properties 
effectively as is done in computer experiments, perturbation 
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theories and variational methods. Reference [15] provides a 
review of these methods. 

(ii) Anisotropic Forces 

The spherically symmetric rn-6-8 potential cannot "be, 
strictly speaking applicable for polyatomic molecules as it 
does not include any angle dependent interations. But in 
predicting phase equilibrium properties it has been shown by 
G-ubbins and Gray [19] that the anisotropic part of the poten- 
tial plays a crucial role. Rear the critical point the cal- 
culated results prove to be very sensitive to the angle- 
dependent terms. In fact in systems that have a liquid- 
liquid -gas line that ends at an upper critical end point, 
which have been studied in detail by Gubbins [19] et al, it is 
found that anisotropic forces are perhaps the most decisive. 
Hence the effect of adding specific angle dependent 
interactions to the m-6-8 rnodeal may be expected to be sig- 
nificant. We carry out a preliminary study of these effects 
in CO2 in this thesis. 

2.6 Anisotropic Potentials 

The earlier works involving calculation of phase 
equilibrium properties used angle dependent potentials 
such as rigid ellipsoids of revolution, sphero-cylindrical 
molecular model (Kihara), rigid spheres containing a point 
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dipole (Keesom) and Stockmayer potential. The details 
of these potentials are given by Hirschf elder et al [20]. But 
these models consider only shape dependent terms, are not 
complete and cannot predict the properties correctly. 

For polarizable, quadrupolar molecules such 
as C 02 » N 2 and F 2 » we mast consider the following general 
potential. 

0(r) = 0(spherical) + 0(quadrupole) + 0(induced dipole) 

+ 0(shape) + 0(anisotropy ) 

The first term is the spherical contribution - the m-6-8 
in our case, the second term represents the electrostatic 
interactions of the permanent quadrupole moments and the third 
term represents the induced-dipole interactions caused by the 
induction effect of the quadrupolar moments and when applicable, 
by the presence of an external electric field. Finally, the 
last two terms depict the anisotropy in the repulsive forces. 

The last two terms may be eneglected for the following 
reasons : 

(i) It is currently impossible to determine independently the 
parameters for the shape part of the potential. 

(ii) The anisotropic contributions have not been worked out in 
enough detail for m~6-8 potential. 
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(iii) The spherical potential does reasonably well and 
hopefully, therefore, incorporates part of these corrections 
into the spherical parameters, especially that of the shape. 

With these assumptions the model potential is modified as ; 
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0 m -6_8 ( M ) + 0 quad (r,Wi,w 2 ) + 0 lnd (r,w 1 ,w 2 ) (2.6.1) 


where w i denotes the set of angles describing the orientation 
of the molecule i with respect to a space-fixed coordinate 
system and r is a vector connecting the centres of mass of 
the two interacting molecules. The spherical term 0 m _g _ 8 ( ! r ! ) 
is given by equation (2.1.2) with |r j = r. 

It is always convenient for linear molecules to write out 
the anisotropic contribution to the potential in terms of the 
relative orientation of the two molecules given by (0-^,©2»0) 

(see Fig. 1) . 


r = distance between the centres of molecules 1 and 2. 

© 1 # ©2 = relative angles made by molecules 1 and 2 with r. 

0 = 01 - 02 = azimuthal angle of molecule 1 with respect to the 
plane containing r and molecule 2. In terns of r, ©^, ©g and 
0 the last two terms in equation (6) become 


0 


quad = ^5 C^*~5 cos 2 ©^^ - 5 cos 2 © 2 - 15 cos'"©^ cos 2 ©2 + 


4r ■ 


2(4cos©-^cos©2 - sin©-^sin ©2 cos0) 


2 n 


( 2 . 6 . 2 ) 
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and 

2 

^ind = ~~~S — (sin^©^ + sin^©2 + 4cos^©^ + 4cos^©2) (2,6.3) 

where a is the mean polarizability and © is the scalar quadrupole 
moment . 

Table 2 provides a and © for an< ^ are ^ken 

from references 2 and 3. 


Table 2 

Gas 

10 26 ©(esu) 

10 (cm^) 

Oxygen 

- 0.4 + 0.1 

1.568 

Nitrogen 

- 1.4 + 0.1 

1.730 

Carbondi oxide 

- 4.3 + 0.2 

2.925 
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CHAPTER 3 

VARIA.TI0MA1 METHOD 


3.1 Introduction 

The essential idea of the variational techniques is to 
expand the thermodynamic properties of the system about a simple 
reference system with a known relationship between its macro- 
scopic and microscopic properties. 

This idea is originally due to Zwanzig [ 21 ]. Using this 
method many works [22, 23] had been carried out based upon a 
simple system of hardspheres. 

Ali Mans o or i and Canfield improved the variational technique 
to calculate the equilibrium thermodynamic properties of simple 
fluids [24] and fluid mixtures [25]; they also extended it to 
apply to the melting process [26]. 

Their variational technique is based on an inequality 
pertaining to the Helmholtz free energy which is discussed below. 
Helmholtz free energy of the original system is calculated by 
variation around the free energy of the reference system based 
on the this inequality and the other thermodynamic properties 
are calculated from the free energy. 
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3.2 Inequality for the Helmholtz Free Energy 

Let us consider a variable x which has a probability- 
distribution function p(x). Then the average of any function 
is defined as 


< f(x)> = f p(x) f(x) dx 

From Taylor series expansion of < f(x) > around f(< x> ) 

2 

< f ( x ) > = f(<x>) + < - ^ - r . < A2_L.il f«< (< x >) + .... 

2 i 

The linear term in x - < x > is absent since < x ~ < x>> = 0. 

The above series gives the two following- inequalities [27] 


(i) < f(x)> * f(< x> ) (3.2.1) 

if f" (x) > 0. 

< ( < > \ 2 

(ii) f (< x> ) f ( <x > ) + — — — ^— > f " ( < x> ) + 

2 '. 


if f 



3 l 


( x ) b 0 


( 3 . 2 . 2 ) 


Let U be the potential energy of the original system and 
U Q be that of the reference system at the same temperature 
and density. Then at equilibrium, we have 

f/f 0 = < exp [-(3(U - U 0 )]> Q 


(3.2.3) 
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where f Q and f^ are the partition functions of the reference 
and the original system respectively. 

Equation (3.2.3) means that f/ f Q is equal to the expe- 
ctation value of exp(-p (U-U 0 ) ) over the probability distribu- 
tion function of the reference system, GuCjp ("PU 0 )/z 0 where 
z Q is the configurational integral of the reference system. 

Substituting x = -p(U - U Q ) (3.2.1) gives 

f/f Q exp(-p <u - U> Q ) (3.2.4) 

and (3.2.2) gives 

f/f Q > exp (-(3 < U - U> 0 ) { 1 + (p 2 /2l )y 2 - 

(3.2.5) 

where y 2 = < (U - U 0 - < U - U> 0 ) 2 > Q 
and y^ = < (U - U Q - < U - U> 0 ) 5 > Q 
¥e know that the Helmholtz free energy A is given by 

A = -kT In f (3.2.6) 

Substituting 


A ( k Q+ <U - U Q > 0 


(3.2.7) 


and 


A $ A n + < U - U n > n - kT In (1 + 4- p 2 y 2 - £ p 3 y^) 


( 3 . 2 . 8 ) 


o 



16 


Inequality (3.2.8) contains y ^ and y^ which are functions of 
second and higher order correlation functions of the reference 
system which are exceedingly complicated [27]. Further in a ther- 
modynamic system with IT particles (of the order 10 , the first 

and second terms in the RHS of (3.2.8) are of the order of NkT. 
while the third term is of the order of kl ln(IT). Consequently 
inequality (3.2.8) will be practically the same as inequality 
(3.2.7). Ali Mansoori et al calculated the Helmholtz free 
energy using (3.2.8) with great difficulty and found that the 
result is not improved even slightly compared to that obtained 
from (3.2.7) but led to some errors. [24]. Hence in this work, 
inequality (3.2,7) is used for further calculations. 


3.3 Representation of <U - U Q 'y c in terms of m-6-8 and 


hards phere Potentials 

< U -U 0 > Q in (3.2.7) can be expressed in terms of the 
m-6-8 potential function by considering only two body interactions 
as follows : 


N 

U = l 

i < a = l 


0 . . (r ■ ■ , . ■ ) 
^lO^iO’ io ' 


and 


IF 


N 

I 


i < 3 


= 1 


0 . . (r • • ) 
*013 v 13 
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For any property x, 


< x> = j~j7 .... J (exp(~pU 0 )/Z 0 ) x dr^ dr 2 .... dr^ dr- 2 


df2 


13 


d.t 


; ‘•( 11 - 1)11 


Therefore, 


< U - U Q > O = I //•••• /(exp(-pu 0 )/z 0 )(u - V di^ ar 2 .... 


dr E dfi 12 ® 13 


d ' 


(H-l)iT 


(3.3.1) 


where k, and E are normalization constants and d 


fl 


10 


by 


dft - • = sin©- d©. sin©. d0- • 

J* . l h-.i 


0 ^10 


xs given 


(3.3.2) 


Z Q , the configurational integral of the reference system is 
given by 


Z 


/ .... J exp(-|3U 0 ) dr-^ dr 2 dr^. 


Separating the angle dependent terns in (3.3.1) 


E 


< U ~ U o > o ~ e/**** / Z q ^ I ^o^ r ij ) dr i dr 2 *** 


dr^ x 


E 

. J . =1 Wij - ^o±j > *12 dfi i 3 •••• dn C3-5 ' 3) 


(E-l)E 


How the angle dependent part can be integrated over all angles 
to give 
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N 


< u - V o = hi ! ■■■ I z~ ex P (-P . < \ , 0o( r ij)) 4r 

1 o i < 0=1 J 


1 dr 2 


... dr^} { 


I 


* (r 13 ) - 8 n 0 o (r i3 ) } 

-L J 


(3.3.4) 


where the integral of the ra-6-8 potential along with angle 

dependent terms (as in equation (3.3.4) ) over d/\ - . is defined 

3-3 

as (r. .). Since 0 (r. -) is independent of t ^ , integration 

X J U X J 

with respect to Q . . will give 8 it 0 (r . ,). Note that 

X J Q X J 

K^_ = 8it and is different from K. 

The radial distribution function g° for the reference 
system is given in terms of second order correlation function as 


g°(r) 


f 2 

7 s 

M(H-l ) 

2 „ J 

P J 


.. jexp(-pn o ) dr- dr^_ ... dr^ 


(3.3.5) 


Representing the integral (3.3.4) by relative vector r^ 2 e ^ c . , 


< U - U > = 

o o 8it 


J fjr— ( exp(“pU 0 ) ] [ i{; (r-^ 2 ) 

8it0 o (r 12 )] dr-j^ dr 2 .... dr^ (3.3.6) 


Combining (3.3.5) and (3.3.6) 
2 


" - u o o = 155 


l «o( r 12> I> ( r 12^ " 8"0 o ( r 12b<ii- 1 dr 2 


V 


r 2 

VT 


^ / g o (*0 r2 [^ (r) - 8it0 (r)] 


cl: 


(3.3.7) 
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where 


l> M ~ 87t0 Q (r)] =J [0(r,Q ) 

ft 

and V is the volume. 


0 o (r)j 


(3.3.8) 


Combining (3.3.7) and (3.3.8) 
r 2 « 


< u - u > 

o o 


|V 7 1,2 dr / L'0( r » ft ) - 0 o ( r )] d ft 
0 ^ 


71 

J 


w 2 00 d n 

= JT / r g (r) / sin« 1 dQ- 
0 

2m 

/ d0 [0(r, ft ) - 0 ] 
o 


sinGg d ©2 


(3.3.9) 


Taking the last integral from (3.3.9) 
2m 

/ [0(r , ) - 0 (r) ] d0 


0 


2m 

1 

0 


/ { [0 sp (r) - 0 Q (r ) ] + ^ M K 

4r 


(1-5 m - 5 cos 2 © 2 - 


2 2 

- 15 cos © p cos ©g + 2(4cos© p cos © 2 


- sin© p sin ©2 cos0) 2 ] - 2 ^£^— (sin^©^ + sin^^) 

8r 

+ 4 cos^© p + 4 cos^Gg] } d0 


2m(0 sp — 0 O ) + ^^*~(l— 5 cos^Gl — 5cos‘ L ©2 — IScos^©^ 

2„ \ 

cos © 2 ) 
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+ — g - 6 - cos 2 ©-, cos 2 ©, + sin 2 ©-, sin 2 ©, 

r 5 1 2 2r 5 12 


Q/vQ ^ ,4 4 4 / 

~ — g- (sin ©^ + sin ©g + 4 cos © 1 + 4cos^'©2) 
4 r 


New, 


TC 

/ 

o 


2tu 

8 in © 2 dOg j 


d0(0(r, jj ) 0 (r) ) 


4Tc(0 a7 , - 0 Q ) + ” 3 “ (2-10 cos^ - -?p - 10 cos 2 © 1 ) 


sp 


+ 2258? cos 2 0. + 30? Sin 2 0, - 2™*! (2 ain 4 

r 9 x r 9 x 4r° ± 


+ 4* 8 cos 4 © 2 . + 


and finally 


71 


% 


2 % 


/ sin© 2 d ©2 J sin© 2 d© 2 J d0[0r, n ) - 0 Q (r)] 


o 


o 


= 8 % - * 0 > 


24naQ‘ 


(3.3.10) 


Equation (3.3.10) contains no quadrupole interaction term since 
it vanishes when averaged over all the orientations. But in 
real systems of CO 2 , the effect of quadrupole interaction is not 
negligible and hence this seems to be a weakness of this 
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approach. The dipole moment of OO2 is zero and hence no 
dipole interaction term is included. Even if a system has 
dipole interaction, the unweighted average of it over all 
the orientations will also vanish as in the case of quadrupole 
int era ct ions . 


This is because in this method only the unweighted 
average of the potential enters the in equality. Physically, 
however, the molecules will tend to spend more time in those 
orientations which give the lowest energy. The average poten- 
tial in these cases may be taken as 




quad 


^gauad [ ex h-^q u aa/"U ] 


;/ 


dw 1 dw 2 


[oxp - (0 quad /kT)] q W;l dw 2 


(3.3.11) 


(3.3.11) cannot be easily evaluated analytically and hence 
(3.3.10) is taken for further calculations. 

Combining (3.3.9) and (3.3.10) 

< u - V 0 > o = ^ f g 0 V)r 2 L(0 sp - 0 O ) - 

0 r (3.3.11) 

Substituting (3.3.11) in (3.2.7) we have 

A < A q + 2jfE£ J g Q (r) r 2 [(0 sp -0 O ) - 2S§^] dr (3.3.12) 


Defining the nondimensional A as 
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A = A/UkT 

(3.3.12) gives 


A * A o + l W sp - 00 - 75“) e 0 M r 2 dr (3.3.13) 

Considering the hard sphere reference system, 



O 

11 

o 


for r £ 1 




CO 


for r < 1 



and 

O 

ii 

to 

o 

for 

r < 1 



where 1 is the hard sphere diameter. Using 

of hard sphere system, (3.3.13) becomes 

the above 

equations 

A* 

A 0 + 2% p 

CO 

p / 

1 

<0 S p - ^e 0 (r) r 2 

dr 

(3.3.14) 


Nondimens ionalising r by defining 
x = r/l 

(3.3.14) becomes 


00 2 

A* A A* + 2n pi 3 / (0 sp - 2 s !-) g 0 (x) x 2 dx (3.3.15) 

1 P 


let us introduce the variational parameter ' c 1 defined as 
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c = l/cr 

0 (m-6-8 potential) as a function of x with the parameter c 

is given by 

0 sp (x) = e[c m (l/cx) m - c 6 (l/cx) 6 - (1/cx) 8 ] (3.3.16) 

where 

c - 6 + 2 v r! m 
c m “ m - 6 a 


„ m - v (m-8) 
°6 = — ~ m 1“ 6 " 


and 

c 8 = V d ® 


With this equation (3.3.15) becomes 


A* 4 : A* + 2tcp p 1? J { [ e(c m (l/cx ) L1 - c g (l/cx) 6 - c 8 (l/cx S ]> 

1 


3a0 2 
7 75 

( cxa) 


g 0 (x) x d: 


(3.3.17) 


Let us define 

2 

V(x) = [c (l/cx) 1 * - c 6 (l/cx) 6 - c 8 (l/'cx)° - ^§-(l/cx) 8 ] 

ea 
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A + 2% p p l?e r V(x) g (x) x 2 dx (3.3.18) 

1 o 

Now the concern is with g (x) for a hard sphere system. Many 
methods are available to evaluate g (x) analytically and nume- 
rical]^ [28] - [32]. Except ¥er the ini' s analytical expression 
[28], all the others are either numerical or in the form of 
complicated integrals which cannot be used conveniently in the 
variational calculations. 

Wertheim [28] has obtained the explicit Laplace transform 
G-(s) of the hard sphere g Q (r) satisfying Percus-Yevick integral 
equation 

00 

G(s) = ^ exp(-sx) g Q (x) x dx (3.3.19) 

1 

g Q (r) obtained from (3.3.19) has two main defects [36] when 
compared with Monte-Carlo on molecular dynamic simulation 
results. 

(i) It predicts a low first peak. This can be corrected by 
defining a new value for a. But in this work as cr is included 
in the variational parameter ’ c* , this correction is not 
necessary. 

(ii) It does, not predict the consecutive peaks and oscillations 
properly. To correct this an additional increment dg(x) should 
be added. We assume that this error does not significantly 
affect the calculation of the thermodynamic properties. 
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Hence (3.3.19) is used in this work without any addition 
of extra terms. 

Defining the nondimensional temperature 

T"' r = 1 /pe 
(3.3.18) becomes 


A* « A* + / 7(x) g (x) x 2 dx (3.3.20) 

u 'T 1 u 

Let us consider V(s) as the Iaplaco inverse of x V(x); 
then — 

CO _ - • 

...-ux'V(xj = * exp(-sx) V(s) ds (3.3.21) 

o 

Substituting (3.3.21) in (3.3.19) and arranging the integrals. 


x V(x) g (x) x dx = / 7(d) Q-(s) d 


as 


o 


(3.3.20) and (3.3.22) give 


A* 1 A* + JV(s) G(s) d£ 


o 


T 


o 


V(s) = l” 1 [x V ( x ) ] 


(3.3.22) 


(3.3.23) 


El- 2 


c - (1/c)n - °6 (1/c)6 fr 


c 8 (l/c) 8 si 

8 sr 


(l/c ) 8 S 6 


ecr 


w 


ST 


(3.3.24) 
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and. 

G-(s) = L[x g Q (x)] 

From Wertheim [26] 

&(s) = a I(s)/12ti[L(s) + S(s) e S j (3.3.25) 

where 

L(s) = 12r)[ (1 + £n)s + (1 + 2rj ) ] (3.3.26) 

S(s) = (1 - n) 2 s 3 + 6 t)(1-t])s 2 + 18ri 2 s - 12 t]( 1 + 2p) 

(3.3.27) 

and 

1 ^ -.3 

Tj = 7t P 1 

1 3 3 

= g 1 it p cr or 

1 * 3 

= TC P C 

* 3 

where p = p cr" is the nondimensional number density. Now 
(3.3.23) can he written as 

00 

A* ^ A* + / V(s) G(s) ds (3.3.28) 

O m 

O 

For a single component thermodynamic syst on 9 a ccording to 
phase rule, if two of the intensive properties of the system 
are given, the other properties will be defined. Of course, 
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along the phase transition line one intensive thermodynamic 
variable of the system will be enough to fix the other variables. 
Then in general if we choose the two independent variables , 
the density and T, the temperature, one could write, 

A = A(p , T) 


and the inequality (3.3.28) becomes 


A* ( 


-X\ 

I ) ' = 


A*( 


* 

P f 


T 


* 


-x- 


c) + F( p , c) 
f 


where 

F( p *, c) = 12r) J V(s) G(s) ds 

o 


(3.3.29) 


The RES of (3.3.29) is a function ofp", T and c whereas the 
LHS is only a f miction of p* and T . This means, that, 
to bring the inequality (3.3.22) closer to equality, RES should 
be required to satisfy the following conditions : 



3 c 


, 1 
T 


8_F 
3 o 


0 



2 

1 3 F > 

— >r~ ' 
T 3 c 


0 


(3.3.30) 


(3.3.31) 


If these conditions are satisfied, there will be a relative 
minimum for the RES of the inequality (3.3.29) with respect to 
c and (3.3.29) becomes 
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where c ^ is the value of e which satisfies the condition 
(3.3.30) and (3.3.31). 
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CHAPTER 4 

CA1CUIATI0HS FOR C0 2 

4.1 Properties of the Hard sphere System 

For hard sphere system, there are several equation of 
state available [34]. The Percus-Yevick equation [35, 36] 
is very simple and is in good agreement with machine 
calculated data. According to P-Y equation 

|~7p = (1 + p + p 2 - |p 5 )/(l - p) 5 (4.1.1) 

(4.1.1) will give the Helmholtz gree energy from the thermody 
namic relation 

EkT ~ T ^-) ^ — 1 (4.1.2) 

o 

and from (4.1.1) and (4.1.2) 

2" i n (i~p) + ii^j + "T ln P/ * 

(4.1.3) 

which will give the free energy of the reference system. 

4.2 Calculations for C0 2 

For spherically symmetric m-6-8 model, the following 
parameters are selected for C0 2 (from table 1). 
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m = 14, v = 1.0, a - 3.68 £, e/lc = 282 °E and d = 1.048729 

Substituting these values in (3.3.24) with proper units, 

V(s) = 4.06395 x 10” 9 s 12 c~ 14 - 0.00206393s 6 c” 8 - 0. 05 54 33 s 4 c~ 6 

( 4 . 2 . 1 ) 

(4.2.1) is true only if interactions are excluded. 

Hanley [7] found the following values are the best for CO 2 
systems when interactions are also included. 

m = 11, v - 0.5, a = 3.775 e/k = 217 °K and d = 1.1259862 

In the following, since we want to include the effect of 
interactions the above values will be used. With this set of 
values (4.2.1) becomes 

V(s) = 1.4230781 x 10“ 5 s 9 c” 11 - 1.976783 x 10" 3 s 6 c~ 8 

- 0.1613395 s 4 c“ 6 (4.2.3) 

To simplify the calculations, let us substitute 

A = V(s) (from (4.2.3)) 

B = [(1 + JrT]) s + (1 + 2rj ) ] 6p 

and 

D = [ (1 + 4-n) s + (1 + 2rj ) ] 12-n + [1 - p) 2 s 3 
+ 6p (l-p)s 2 + 18p 2 s - 12p (l+2p)]e S 
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Now (3.3.28) becomes 


A*< a*+ — 

A * A Q + p / D d ° 


(4.2.4) 


A* is given by (4.1.3). Differentiating A* with, respect to c 


8 A o _ 9 A o 3 jq, 

9 0 9 tj 9 c 

(2 - n - } N 2 ) * 2 

3 — (« P G ) 

(1 - nr 

The condition (3.3.30) becomes 


„ 1 2 

I---?-- .£L U * 0 2) + 2_ 7 1 fD ' 4 SB + B CLA > AB dD-, ds = 0 

(1-ri)^ V P ^ T* 4> D^ L dC dC * d ° J 

(4.2.5) 


where 

dA _ 
dC ~ ‘ 


■15.653859 x 10~ 5 s 9 c” 12 + 15.814264 x 10“ 3 s 6 c“ 9 


+ 0.968037 s 4 c 7 


( 4 . 2 . 6 ) 


|§ = 3ti p y c 2 (i+s+4n+5n) 


(4.2.7) 


= n;p* o 2 (6+6s~6e S + 3s 2 e s - s 5 e s ) + 6%r\ p* c 2 (4+s-4e S + 


„ s 2 s , s^ Sx 
+ 3 s e -so + r~ e ) 


( 4 . 2 . 8 ) 


The integral in (4.2.5) cannot be evaluated analytically after 
substituting (4.2.6), (4.2.7) and (4.2.8) and hence it was eva- 
luated numerically. 


5 .1. i i 

ENIh'AL. I 


JSfe, £% 



THE CUBE OF ' 
PARAMETER 'c ' 
AGAINST R E D U i 
FOR VARIOUS R 

perature 



REDUCED DENSITY * * 


Variational parameter 
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( i ) The Helmholtz free energy 

The inequality (4.2.4) will give the Helmholtz free energy 

when the integral of RHS is minimized; this integral has to be 

evaluated numerically. Similarly the integral in (4.2.5) has to 

be evaluated numerically to find the c which, satisfies this 

equation. That particular c (which is c . ) should be substituted 

nun/ 

into (4.2.4) to get the 'correct' Helmholz free energy for a 
particular temperature and density condition. 

It is very difficult to check the condition (3.3.31) which 

involves the second derivative of A* with respect to ' c' . But 

it can be verified indirectly. The ’ c' which satisfies the 

condition (3.3.30), (equation (4.2.5), may either give maximum or 

minimum value of A* , can be evaluated from (4.2.4) and (4.2.5). 

From the trend of the curve 'A* * versus ' c ! in the neighbourhood 

of c . . one can immediately conclude whether c . gives a 
min/ 17 mm D 

maximum or minimum. 

By evaluating the values of variational parameter ' c' which 
give the minimum value of A* , it is evident that there is only 
one such 1 c' for every T* + . 1 11 otherwords, RHS of 

(4.2.4) has only one minimum with respect to ' c’ . 

The effect of temperature is studied from the calculation. 

It Is found that at high temperatures, the variational parameter 
which is the ratio of hard sphere cut off diameter 'A* to 1 cr* Is 
nearly unity (0.93 to 0.99) and it is also found that at there 
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tempera tur e the ratio is almost independent of density. 

But at very low temperatures, well Below the critical 
tempera tured the ratio is found to he a strong function of 
density. Also the deviation from unity is more compared to 
the high temperature and the values are found to he greater 
than unity. 

This shows that there is great deviation from the hard- 
sphere behaviour at very high temperatures because, when low 
energy collisions occur, the interacting system is strongly 
subject to the potential well. This fact is well established 
byHHirsch fields through quantum mechanical calculations. 

The integrals in all these equations have limits zero 
to infinity and hence xfhen they are evaluated, a convergence 
criterion should be considered for each thermodynamic condition 
to get the upper limit. The lower limit also gives computational 
troubles. When the variable s takes the values near zero, either 
under or overflow occurs. 

In this work, these difficulties were minimized by using 
Gauss -la guerre integration method by fitting the curve to a 
polynomial having 14 zeros. This method seems to be more 
accurate than Simpson' s method with reasonably small increment 
in s. Further, computational time is found to be less than 
l/lOOth of the required for Simpson’s method when the increment 
used in the latter method is such that it gives the same 
result as the former method. 
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Using these values of c a±ri values of A* are obtained 
for various sots of nondimens ional temperature (I*' ) and 

X* 

nondiniensional density ( ? ) . 


(ii) Equation of State 

(4.2.9) 

¥e know that 

§V^T,N = (4.2.10) 

and 

f= U/V (4.2.11) 


Therefore 


^9V^t,u 


f 2 /u 


( 4 . 2 . 12 ) 


Substituting (4.2.12) ain (4.2.10) and combining (4.2.9) 

2 

= P 


A 


_) -P- 

p 'T,U U 


Multiplying both sides by V/hkT, 

f _ PV „ 

^ ; T,N ~ NkT “ 

As p*= p a 5 , (4.2.13) becomes 




(4.2.13) 


Z 


(4.2.14) 
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Equation (4.2.14) shows that the product of p* with the 
differential coefficient of A" with respect to P will give 
the compressibility factor. 

¥e have evaluated numerically for various p* and . 

For this purpose, fourth order difference method is used and is 

found to be very accurate. The derivative is evaluated conti- 

?«■ 

nuously as a function of f with a regular interval and also 
at discrete points. The second method is less time consuming 
and the accuracy is adjustable without any loss of computer time. 
Comprenibility isotherms have been evaluated for wide range 
of temperatures and graphs are drawn for some selected tempera- 
tures and are shown in Fig. 3 and Fig. 4 and also are given in 
Ta ble 3 • • 

The presures calculated from this equation of state are 
compared with experimental result in Fig. 5 and the agreement is 
found to be good. At high temperatures the error is well within 
5 percent for practical pressure range. However, at very high 
pressure, the error even exceeds 5 percent. But, it is to be 
mentioned that at all densities and temperatures, the prediction 
of this method are more accurate [38, 39] than the previous 
methods using I*J model. The spherical m-6-8 potential calcula- 
tions are also compared with these results in Fig. 5 and it is 
found that the results obtained with the inclusion of dis- 
persion on terms are better at 375°C. 
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(iii) Phase equilibrium calculations 


When two phases are at equilibrium, the temperatures, the 
pressures and the chemical potentials of the two phases should 
be equal. Considering the liquid and vapour phases, we have 

P L = P V (4.2.15) 
(A+PV) l = (A+PV) V (4.2.16) 


Defining the nondimensional pressure as 


Per 

e 


(4.2.17) 


Condition (4.2.15) becomes 


■X* __X- 

P L = P V 


(4.2.18) 


Dividing (4.2.16) by WkT, 


(A* + Z), = (A* + Z) 


V 


(4.2.19) 


P can be related to Z as 

P* = Z P * T* (4.2.20) 

For a particular temperature we have A* , Z and hence P* as a 
function of P *. In the system of two equations (4.2.18) and 
(4.2.19), the only unknowns are the two densities of the co- 
existing phases which can be numerically solved using the 
calculated values of A* and P for a particular T* as a 

function p , 
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At low densities, the interval in P* is not sufficiently 
close to get a reasonably accurate solution to the system of 
equations (4.2.18) and (4.2.19) and hence the interval which 
gives accurate solution for higher density range should be 
further divided into subintervals for low density range. As 
our interest is mainly on the critical region, the p ha se 
equilibrium calculations are carried out near the critical 
point and a plot of temperature versus co-existing densities 
of liquid and vapour is made as shown in Pig. 7. Phe logarithm 
of the vapour pressure predicted from this procedure is also 
plotted against inverse P* as shown in Pig. 8. 

Phe phase equilibrium and vapour pressure calculations 
prove that this method is quite satisfactory near the critical 
region. Phe phase equilibrium curves as shown in Fig. 7 are 
very close to the curves using the law of rectilinear diameter 
and the experimental values. While the saturated vapour locus 
is predicted with negligible error, rthe saturated liquid locus 
is predicted with a maximum error of about G percent. 

It may be anticipated that the errors in these calculations 
will be reduced considerably if one takes the quadrupole 
interaction into account by taking the weighted average of this 
interaction over all the orientations. It is suspected that 
quadrupole interactions may play an appreciable part in the 
interaction between OOg molecules [43]. 
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Another improvement can be made in these calculations by 
talcing a better reference system. When the reference system 
approaches the original system, the inequality (4.2.4) will 
tend to an equality. 

(iv) Other thermodynamic properties 

Thermodynamic properties like entropy, enthalpy and 
internal energy can be readily evaluated from A and Z which 
are already evaluated. 

From thermodynamic relations, 

6 a/ 3 T) v ^ w = -s 

*)£ 

From the definition of A, 

A = A* * NkT, 

-g-7£ = '|| NkT + DlcA* = -S (4.2.21) 

Defining S'* = S/Nk, 

-S* = T* f~+ A* (4.2.22) 

Internal energy U can be related to A by, 

U = A + TS (4.2.23) 

Dividing by MIcT on both sides and defining 

U* = A* + S* 


(4.2.24) 
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Combining (4.2.22) and (4.2,24) 

u* = -I* 4 t* (4.2.25) 

and from (4.2.24) 

S* = U* -A* (4.2.26) 

-Jf -)f 

Prom the equations (4.2.25) and (4.2.26) U and S can be 
evaluated. 

The nondimensional enthalpy H* is defined by 

H* = Ii/NlcT (4.2.27) 

and from the thermodynamic relation 

H* = A + PV + TS (4.2.28) 

¥e got 

H* = A* + Z + S* (4.2.29) 

Equation (4.2.29) will give one H* for each value of * and T* 
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CHAPTER 5 

conclusion 

The pressures calculated from this equation of state are 
compared with experimental results (Rig. 3,4 and 5) and the 
agreement is found to be good. At high temperatures there is 
negligible error and at low temperatures the error is well 
within 5 percent for practical pressure range. 

The phase equilibrium and the vapour pressure calculations 
(Pig. 7 and 8) prove that this method is excellent near critical 
region. While the saturated vapour locus is predicted without 
any error, the saturated liquid locus is predicted with a maxi- 
mum error of 6 percent. 

The error in the liquid locus and in the values of pressures, 
at low temperatures lead to two conclusions t (i) The qudrupole 
interaction term becomes zero in our calculations because only the 
unweighted average of the quadrupole moment term is considered. 
This has been already discussed in the calculation part. The 
absence of the quadrupole interaction term in the equation of 
state affects these calculations much at low temperatures and at 
higher pressures. It is suspected that the quadrupole interaction 
may play an appreciable in the interaction between OO2 molecules 
at low temperature and high pressure [41], (ii) Ihe simple and 
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angle dependent hard sphere potential is taken as the reference 
system but on the otherhand the original system is angle dependent 
and containing more terms. Owing to this inadequency of the 
reference system, we can not expect the inequality for Helmholtz 
free energy, from which all the calculations are made, will ■ 
tend to an equality even at the minimum value of A 

The spherical m-6-8 potential calculations with the best 
value available for m are compared with this calculations. It 
is also compared with L.J. potential calculations. It is 
concluded that this work leads to bettor results than the other 
two at all the thermodynamic conditions. Hence we can conclude 
that the addition of angle dependent terms improves the accuracy 
in the calculation of thermodynamic properties and gives satisfa- 
ctory results. 
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Table 3 

Relation between dimensional and nondimensional PVT 


Property Relation to non- Relation for 

dimensional property carbondi oxide 


1. Temperature in T 

°K 

2. Pressure in P 

atmospheres 

3. Density in gm/co P 




M 


N <? 


T = 217 T* 


P = 549.378 P* 

P = 1.3579661 p* 
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PVT data for CC>2 

T * = 0.1 


DENSITY p * 

PRESSURE P 

DENSITY p * 

PRESSURE P 

0.02 

- 0.00045 

0.52 

- 2.01660 

0,04 

- 0.00545 

0.54 

- 2.18110 

0.06 

- 0.01551 

0.56 

- 2.34990 

0.08 

- 0.03086 

0.58 

- 2.52232 

0.10 

- 0.05173 

0.60 

- 2.69764 

0.12 

- 0.07835 

0.62 

- 2.87492 

0.14 

- 0.11093 

0.64 

- 3.05320 

0.16 

- 0.14969 

0.66 

- 3.23141 

0.18 

- 0.19480 

0.68 

- 3.40819 

0.20 

— 0. 24646 

0.70 

- 3.58231 

0.22 

- 0.30482 

0.72 

- 3.75207 

0.24 

- 0.37003 

0.74 

- 3.91584 

0.26 

- 0.44220 

0.76 

- 4.07163 

0.28 

- 0.52144 

0.78 

- 4.21738 

0.30 

- 0.60783 

0.80 

- 4.35076 

0.32 

- 0.70140 

0.82 

- 4.46905 

0.34 

- 0.80217 

0 . 84 

- 4.56957 

0.36 

- 0.91014 

0.86 

- 4.64910 

0.38 

- 1.02523 

0.88 

- 4.70409 

0.40 

- 1.14737 

o.so 

- 4.73091 

0.42 

- 1.27640 

0.92 

- 4.72521 

0.44 

- 1.41215 

0.94 

- 4.68266 

0.46 

- 1.55438 

0.96 

- 4.59770 

0.48 

- 1.70277 

0.98 

- 4.46546 

0.50 

- 1.85699 

1.00 

- 4.27946 
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PVT data for CCu 

* ^ 
T = 0.2 


DENSITY p * 

PRESSURE P 

0.02 

- 0.00148 

0.04 

- 0.00118 

0.06 

- 0.00874 

0.08 

- 0.02137 

0.10 

- 0.03934 

0 » 12 

- 0.06283 

0.14 

- 0.09200 

0 . 16 . 

- 0.12705 

0.18 

- 0.16815 

0 . 20 

- 0.21546 

0.22 

- 0.26911 

0 . 24 

- 0.32922 

0.26 

- 0.39587 

0 . 28 

- 0.46914 

0.30 

- 0.54905 

0.32 

- 0.63562 

0.34 

- 0.72882 

0.36 

- 0.82857 

0.38 

- 0.93479 

0.40 

- 1.04732 

0.42 

- 1.16596 

0.44 

- 1. 29044 

0.46 

- 1.42049 

0.48 

- 1.55569 

0.50 

- 1.69562 


DENSITY p * 

PRESSURE P *- 

0.52 

- 1.83975 

0.54 

- 1.98748 

0.56 

- 2.13811 

0.58 

- 2.29082 

0.60 

- 2.44477 

0.62 

- 2.59886 

0 . 64 

- 2.75197 

0.66 

- 2.90284 

0.68 

- 3.04993 

0.70 

- 3.19175 

0.72 

- 3.32638 

0.74 

- 3.45192 

0.76 

- 3.56604 

0.78 

- 3.66637 

0.80 

- 3.75022 

0.82 

- 3.81438 

0.84 

- 3.85578 

0.86 

- 3.87060 

0.88 

- 3.85474 

0.90 

- 3.80387 

0.92 

- 3.71294 

0.94 

- 3.57667 

0.96 

- 3.38864 

0.98 

- 3.14282 

1.00 

- 2.83146 
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T * = 0.5 


DENSITY p * 

PRESSURE P * 

DENSITY P * 

PRESSURE P * 

0.02 

0.05402 

0.52 

- 1.765608 

0.04 

0.005072 

0.54 

- 1,909818 

0.06 

- 0.00198 

0.56 

- 2.05652 

0.08 

- 0.011904 

0.58 

- 2.204058 

0.10 

- 0.02694 

0.60 

- 2.5517 

0,12 

- 0.047196 

0.62 

- 2.498166 

0.14 

- 0.072954 

0.64 

- 2.641728 

0.16 

- 0.1044 

0.66 

- 2.780712 

0.18 

- 0.141696 

0.68 

- 2.912712 

0.20 

- 0.1851 

0.70 

- 5.0566 

0.22 

- 0.254762 

0.72 

- 5.148652 

0.24 

- 0.290756 

0.74 

- 5.24675 

0.26 

- 0.555418 

0.76 

- 5.527204 

0.28 

- 0.422772 

0.78 

- 5.586682 

0.50 

- 0.49896 

0.80 

- 5.42096 

0.52 

- 0.582048 

0.82 

- 5.425504 

0.54 

- 0.67218 

0.84 

- 5.594692 

0.56 

- 0.769176 

0.86 

- 5.525298 

0.58 

- 0.875554 

0.83 

- 5.204696 

0.40 

- 0.98456 

0.90 

- 5.05102 

0.42 

- 1.101114 

0.92 

- 2.794224 

0.44 

- 1.226808 

0.94 

- 2.484984 

0.46 

- 1.554608 

0.96 

- 2.08944 

0.48 

- 1.487576 

0.98 

- 1.58615 

0.50 

- 1.6245 

1.00 

- 1.0415 



47 


T * = 0.4 

DENSITY P * 

PRESSURE P * 

DENSITY P * 

PRESSURE P * ~ 

0.02 

0.005336 

0.52 

- 1.560208 

0.04 

0.007328 

0.54 

- 1.684368 

0.06 

0.004776 

0.56 

- 1.81216 

0.08 

- 0.002464 

0.58 

- 1.935576 

0.10 

- 0.01456 

0.60 

- 2.05656 

0.12 

- 0.031632 

0.62 

- 2.17372 

0.14 

- 0.053928 

0.64 

- 2.285312 

0.16 

- 0.081536 

0.66 

- 2.388936 

0.18 

- 0.114624 

0.68 

- 2.482816 

0.20 

- 0.15336 

0.70 

- 2.56396 

0.22 

- 0.197912 

0.72 

- 2.62944 

0.24 

- 0.248352 

0.74 

- 2.676136 

0.26 

- 0.304824 

0.76 

- 2.700128 

0.30 

- 0.36736 

0.78 

- 2.69724 

0.32 

- 0.51072 

0.80 

- 2.66272 

0.34 

- 0.5916 

0.82 

- 2.587592 

0.36 

- 0.678672 

0.84 

- 2.483792 

0.38 

- 0.7714 

0.86 

- 2.305832 

0.40 

- 0.86992 

0.88 

- 2.097216 

0.42 

- 0.974232 

0.90 

- 1.83384 

0.44 

- 1.082752 

0.92 

- 1.500704 

0.46 

- 1.197288 

0.94 

- 1.108072 

0.48 

- 1.314432 

0.96 

- 0.638976 

0.50 

- 1.4362 

0.98 

- 0.08624 



1.00 

+ 0.554 



48 


T * = 0.5 


DENSITY p 

PRESSURE P * 

DENSITY p * 

' PRESSURE P 

0.02 

0.00726 

0.52 

- 1.35226 

0.04 

0.01158 

0.54 

- 1.45611 

0.06 

0.01152 

0.56 

- 1.5582 

0.08 

0.00696 

0.58 

- 1.65909 

0.10 

- 0.0022 

0.60 

- 1.7535 

0.12 

- 0.01614 

0.62 

- 1.785 

0.14 

- 0.03493 

0.64 

- 1.92352 

0.16 

- 0.058732 

0.66 

- 1.99584 

0.18 

- 0.08775 

0.68 

- 2.05326 

0.20 

- 0.1219 

0.70 

- 2.0965 

0.22 

- 0.16148 

0.72 

- 2.1186 

0.24 

- 0.20652 

0.74 

- 2.12084 

0.26 

- 0.25615 

0.76 

- 2.09722 

0.28 

- 0.31276 

0.78 

- 2.0455 

0.30 

- 0.37425 

0.80 

- 1.9564 

0.32 

- 0.4408 

0.82 

- 1.83516 

0.34 

- 0.51323 

0.84 

- 1.66782 

0.36 

- 0.58914 

0.86 

- 1.45641 

0.38 

- 0.61279 

0.88 

- 1.18448 

0.40 

- 0.7592 

0.90 

- 0.86805 

0.42 

- 0.85071 

©.92 

- 0.47564 



0.94 

- 0.02162 

0.44 

- 0.94572 



0.46 

- 0.10442 

0.96 

0.5188 



0.98 

1.12945 

0.48 

- 1.14504 



0.50 

- 1.24825 

1.00 

1.839 


49 


* 

I =0.6 


DENSITY P * 

PRESSURE 

0.02 

0.00919 

0.04 

0.01584 

0.06 

0.01824 

0.08 

0.01634 

0.10 

0.01008 

0.12 

- 0.00069 

0.14 

- 0.01604 

0.16 

- 0.03612 

0.18 

- 0.06096 

0.20 

- 0.09069 

0.22 

- 0.12535 

0.24 

- 0.16494 

0.26 

- 0.20951 

0 . 28 

- 0.25901 

0.30 

- 0.31307 

o 

• 

rv > 

- 0.37250 

0.34 

- 0.43538 

0.36 

- 0.50233 

0.38 

- 0.57304 

0.40 

- 0.64714 

0.42 

- 0.72417 

0.44 

- 0.80363 

0.46 

- 0.88493 

0.48 

- 0.96736 

0.50 

- 1.05013 


T " T 


DENSITY P 

PRESSURE P 

0.52 

- 1.13236 

0.54 

- 1.21300 

0.56 

- 1.29093 

0.58 

- 1.36483 

0.60 

- 1.43328 

0.62 

- 1.49462 

0.64 

- 1.54707 

0.66 

- 1.58858 

0.68 

- 1.61688 

0,70 

- 1.62952 

0.72 

- 1.62364 

0,74 

- 1.59619 

0.76 

- 1.54370 

0.78 

- 1.46237 

O.SO 

- 1.34801 

0.82 

- 1.19595 

0.84 

- 1.00058 

0.86 

- 0.75663 

0.88 

- 0.45733 

0.90 

- 0.09580 

0.92 

+ 0.33623 

0.94 

0.84722 

0.96 

1.46155 

0.98 

2.12607 

1.00 

2.97015 
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* 

T = 0,7 


DENSITY p * 

PRESSURE P * 

DENSITY p * 

PRESSURE P 

0.02 

0.01111 

0.52 

- 0.95550 

0.04 

0.02008 

0.54 

- 1.01939 

0.06 

0.02496 

0.56 

- 1.07913 

0.08 

0.02573 

0.58 

- 1.13333 

0.10 

0.02234 

0.60 

- 1.18041 

0.12 

0.01473 

0.62 

- 1.21857 

0.14 

0.00275 

0.64 

- 1.24583 

0.16 

- 0.01356 

0.66 

- 1.26003 

0.18 

- 0.03440 

0.68 

- 1.25861 

0.20 

- 0.05970 

0.70 

- 1.23897 

0.22 

- 0.08948 

0.72 

- 1.19795 

0.24 

- 0.12388 

0.74 

- 1.13226 

0. 26 

- 0.16262 

0.76 

- 1.03811 

0.28 

- 0.20579 

0.78 

- 0.91136 

0.30 

- 0.25328 

0.80 

- 0.74745 

0.32 

- 0.30485 

0.82 

- 0.54109 

0.34 

- 0.36026 

0.84 

- 0.28678 

0.36 

- 0.41933 

0.86 

0.02704 

0.38 

- 0.48145 

0.88 

0.38423 

0.40 

- 0.54709 

0.90 

0.83317 

0.42 

- 0.61373 

0.92 

1.34806 

0.44 

- 0.68193 

0.94 

1.89765 

0.46 

- 0.75004 

0.96 

2.56773 

0.48 

- 0.82028 

0.98 

3.27639 

0.50 

- 0.88876 

1.00 

4.12241 



51 


T* = 0.8 


DENSITY p * 

PRESSURE P 

DENSITY p 

PRESSURE P * 

0.02 

0.01304 

0.52 

- 0.77865 

0.04 

0.02433 

0.54 

- 0.82577 

0.06 

0.03167 

0.56 

- 0.86734 

0.08 

0.03510 

0.58 

- 0.90183 

0.10 

0.03459 

0.60 

- 0.92755 

0.12 

0.03008 

0.62 

- 0.94251 

0.14 

0.02152 

0.64 

- 0.94461 

0.16 

0.00885 

0.66 

- 0.93146 

0.18 

- 0.00796 

0.68 

- 0.90034 

0.20 

- 0.02886 

0.70 

- 0.84842 

0.22 

- 0.05396 

0.72 

- 0.77103 

0.24 

- 0.08309 

0.74 

- 0.67135 

0.26 

- 0.11625 

0.76 

- 0.52887 

0.28 

- 0.15319 

0.78 

- 0.36433 

0.30 

- 0.19399 

0.80 

- 0.15533 

o 

s 

ro 

- 0.23825 

0.82 

0.09627 

0.34 

- 0.28572 

. 0.84 

0.39488 

0.36 

- 0.33611 

0.86 

0.74577 

0.38 

- 0.38902 

0.88 

1.15446 

0.40 

- 0.44394 

0.90 

1.62679 

0.42 

- 0.50053 

0.92 

2.16274 

0.44 

- 0.55744 

0.94 

2.79446 

0,46 

- 0.61497 

0.96 

3.46770 

00 

e 

o 

- 0.67221 

0.98 

4 . 28143 

0.50 

- 0.72544 

1.00 

5.13329 


52 


T* = 0.9 


DENSITY P* 

PRESSURE P* 

DENSITY p * 

PRESSURE P* 

0.02 

0.01496 

0.52 

-0.59851 

0.04 

0.02856 

0.54 

-0.62990 

0.06 

0.03837 

0.56 

-0.65265 

0.08 

0.04445 

0.58 

-0.66912 

0.10 

0.04680 

0.60 

-0.67434 

0.12 

0.04539 

0.62 

-0.66714 

0.14 

0.04021 

0.64 

-0.64574 

0.16 

0.03119 

0.66 

-0.60746 

0.18 

0,01836 

0.68 

-0.54789 

0.20 

0.00174 

0.70 

-0.47103 

0.22 

-0.01862 

0.72 

-0.36108 

0.24 

-0.04266 

0.74 

-0.23130 

0.26 

-0.07020 

0.76 

-0.05908 

0.28 

-0.10121 

0.78 

0,13980 

0.30 

-0.13537 

0.80 

0.38418 

0.32 

-0.17244 

0.82 

0.67261 

0.34 

-0.21222 

0.84 

1.00986 

0.36 

-0.25406 

0.86 

1.40087 

0.38 

-0.29744 

0 . 88 

1.85088 

0.40 

-0.34294 

0.90 

2.36741 

0.42 

-0.38880 

0.92 

2.94752 

0.44 

-0.43424 

0.94 

3.62095 

0.46 

-0.47889 

0.96 

4.35860 

0.48 

-0.52260 

0.98 

5.19272 

0.50 

-0.56190 

1.00 

6.12930 


53 


* 

T = 1.0 


DENSITY P * 

PRESSURE P * 

DENSITY p * 

PRESSURE P 

0.02 

0.01688 

0.52 

- 0.42354 

0.04 

0.03280 

0.54 

- 0.43810 

0.06 

0.04507 

0.56 

- 0.44425 

0.08 

0.05379 

0.58 

- 0.44214 

0.10 

0.05900 

0.60 

- 0.42753 

0.12 

0.06067 

0.62 

- 0.39982 

0.14 

0.05883 

0.64 

- 0.35593 

0.16 

0.05344 

0.66 

- 0.29320 

0.18 

0.04455 

0.68 

- 0.21199 

0.20 

0.03221 

0.70 

- 0.10497 

0.22 

0.01649 

0.72 

0.02961 

0.24 

- 0.00250 

0.74 

0.19311 

0.26 

- 0.02460 

0.76 

0.38732 

0.28 

- 0.04962 

0.78 

0.62861 

0.30 

- 0.07732 

0.80 

0.90049 

0.32 

- 0.10741 

0.82 

1.22417 

0.34 

- 0.13960 

0.84 

1.59828 

0.36 

- 0.17317 

0.86 

2.02764 

0.38 

- 0.20815 

0.88 

2.51375 

0.40 

- 0.24378 

0 . 90 ' 

3 . 07188 

0.42 

- 0.27864 

0.92 

3.70504 

0.44 

- 0.31333 

0.94 

4.41309 

0.46 

- 0.34602 

0.96 

5.19052 

0.48 

- 0.37603 

0.98 

6.09348 

0.50 

- 0.40185 

1.00 

7.05255 


54 


* 

T 


= 2.0 


PVT DATA POP CO 2 

T* = 3.0 


DENSITY p* PRESSURE P* DENSITY p* PRESSURE P* 


0.10 

0.18 

0.10 

0.30 

0.20 

0.328 

0.20 

0.624 

0.30 

0.480 

0.30 

1.008 

0.40 

0.688 

0.40 

1.56 

0.50 

1.06 

0.50 

2.37 

0.60 

1.776 

0.60 

3.708 

0.70 

3.08 

0.70 

5.796 

0.80 

5.344 

0.80 

9.072 

0.90 

9.072 

0.90 

14.256 

1.00 

14.92 

1.00 

22.2 

1.10 

23.76 

1.10 

34.254 

1.20 

37.152 

1.20 

52.704 

1.30 

58.188 

1.30 

80.652 

1.40 

90.44 

1.40 

123.48 
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